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Abstract 



A scenario of testing quantum mechanics in neutrino oscillation is presented. 
The quantum mechanics violation (QMV) that is motivated by arguments based 
on quantum gravity is investigated in neutrino system. It is found that the evolu- 
tion equation of density matrix including QMV effect is analytically resolvable for 
neutrino propagate in vacuum or in matter under adiabatic approximation. The 
analytical formulas have been derived. Some bounds on the related parameters 
have been obtained from neutrino experiments. 

PACS number(s): 14.60.Pq, 03.65.Bz 



1 Introduction 

It is most essential and important issue that whether neutrino has nonzero mass and 
whether there is neutrino oscillation because which is related to the theoretical struc- 
ture of standard model (SM) of particle physics and is also related to astrophysics and 
cosmology. Recent experiments^-]^ have shown that there are signals of neutrino os- 
cillation. If it is confirmed, neutrino system may not only be a window beyond SM but 
also be a fruitful system for testing quantum mechanics and probing discrete symmetries 
of Nature as in K^K^ system. 

The suggestion that quantum coherence might be lost at the microscopic level was 
made by HawkingQ, which suggested that asymptotic scattering should be described 
in term of a superscattering operator /S*, relating initial and final density matrices, that 
does not factorize as a product of S— and matrix elements 



The loss of quantum coherence was thought to be a consequence of microscopic quantum- 
gravitational fluctuations in the space- time background. Ellis, Hagelin, Nanopoulos and 
Srednicki (EHNS)0 then pointed out that if Eq.(l) is correct for asymptotic scattering, 
there should be a corresponding effect in the quantum mechanics Liouville equation that 
describes the time evolution of the density matrix p(t) 



Pout —^Pi 



(1) 



dtp{t)=i[p,H]+i5 fip 



(2) 



The extra term may evolve a pure state into mixed state and result in quantum mechanics 
violation (QMV). The 5 /H was parameterized for simple two level system by EHNS[0] 



and was used in neutron and K^K^ systems. Then the QMV was through studied by 
the authors of references Il8[] , |P| , [|T0[| for K^K^ system, and by the authors of Ref.|]lll for 



BB system. Some restrictions on QMV parameters have been given. 

A beam of two flavor neutrinos is a simple two-state system. If mass eigenstate is 
not degeneracy with weak interaction eigenstate, then mixing between two flavor will 
exist and neutrino oscillation will occur. As in K^K^ system, the quantum mechanics 
violation originate from quantum gravitational fluctuation may also exist in neutrino 
system. We investigate QMV in neutrino system in term of the formalism proposed by 
EHNS[0] and find that neutrino oscillation probability can be modified by QMV effect. 
Thus the neutrino oscillation experiments may present a precise testing for quantum 
mechanics. 

The main goal of this paper is that establish a framework of QMV in neutrino system, 
derive the detailed formulas of neutrino oscillation probability including QMV effect, 
and present some bounds on QMV parameters. For completeness we will first derive 
the formulas of neutrino oscillation probabilities described by density matrix formalism 
within conventional quantum mechanics in section II. The QMV for neutrino oscillation 
in vacuum and in matter will be investigated in sections III and IV respectively. We will 
apply the formulas derived in previous section to current neutrino experiments and draw 
out certain restrictions on QMV parameters in section V. Our conclusions are presented 
in section VI. 



2 Neutrino oscillation described in term of density- 
matrix formalism 

We now study a system consisted of two- flavor neutrinos in vacuum. The Hamiltonian 
of this system, in mass eignstate, is diagonal, H = diag{Ei, E2), here Ei and E2 are 
energy eigenvalues of neutrinos. In some practical problems, for instance solar neutrino 
detection or atmospheric neutrino observation, rrij <^ Ej are always valid, then Ei ^ 
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E + i = 1,2, are good approximations, here E =\ p \ is the magnitude of three 
momentum of neutrinos. We describe the time evolution of the neutrino system by 
Liouville equation in conventional quantum mechanics 

lp(t) = t[p, H], (3) 

where p is the density matrix of the system. Following EHNS[7], we expand the Hamil- 
tonian and density matrix in Pauli matrix basis 

H = ^h^aa„ p = ^pi3ap, a, (3 = 0,1,2,3, (4) 

where cxo is unit matrix and aj {j = 1,2,3) are Pauli matrices. The evolution of the 
components of desity matrix obey 

dtPa = Kppp. (5) 
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We assume that there are only electron neutrinos in the system at the moment t = 0. 
Considering this initial condition and solving the differential equation of Eq.(5), we get 
the density matrix 

cos'^9 \sin26e'^ 
|sm2^^e"*2E* siri^O 



pit) 



(6) 



where the 9 is the vacuum mixing angle of neutrinos. The probability that electron 
neutrino oscillate to muon neutrino is given by computing the expection value of Tr{Op) 
of an observable 



0{v, 



-\sin2Q 
-^sin26 cos^O 



It is easy to check that the oscillation probability 



(7) 



(8) 



is consistent with that obtained from solving propagation equation of neutrino [|12 
If neutrino propagate in matter, the Hamiltonian is corrected as [ITS 



H = E + 



TTln 



4E 



where 



M2 = - 
2 



~Acos2e + 2A Asin2e 
Asin29 Acos2e 



A = 2V2GFNeE 



(9) 



(10) 



the Ne and Nn are the number densities of electrons and neutrons in matter. The 
Hamiltonian can be diagonalized by an unitary matrix U{6), here the 6 is effective 
mixing angle of neutrinos in matter, which is determined by 



tan2e 



Asin29 



Acos29 - A 

The diagonalized Hamiltonian can be expanded in Pauli matrix basis as 

rid = — o-Q - -Acts, 



here 



S = 2[E - ^GpN^ + ^{ml + ml + A)], 



X = —^{Acos2e - A)2 + A'^smne. 



(11) 

(12) 

(13) 
(14) 



In this situation, the hap in Eq.(5) becomes 



1 ° 
















A 








-A 








I 








; 



(15) 
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In general, it is difficult to solve analytically the evolution equation of the components 
of density matrix, since the A is a complicated function of time. However, it is shown 
(see Appendex A) that if adiabatic condition]!^ is valid, then the approximation 



is available, the term being proportional to (dX/dt) can be neglected during solve the 
differential equation of density matrix. The evolution equation of density matrix can be 
analytically solvable. We still assume that there is only electron neutrino in system and 
the effective mixing angle is 6 at the moment t = 0. Under this initial condition the 
desity matrix can be evaluated out 

. ^ _ / cos'^6 \sin29exp{i Jq Xdt') \ , . 

' ~ \ IsinOexpi-i Jq Xdt') sin'^O ) ' ^ ^ 

If neutrinos are detected in vacuum, the oscillation probability is calculated as 

1 - 1 - /■* 

Pii^e^ I'l,) = -{I- cos29cos2e) - -sin2esm2ecos{ Xdt'). (18) 

2 2 Jo 



Which is just the well known MSW solution[15| of solar neutrino problem. 



3 QMV for neutrino oscillation in vacuum 

We now introduce QMV in neutrino oscillation. Following EHNS0 the extra term in 
Eq.(2) may be parameterized as a symmetrical 4x4 matrix /i^^, then the evolution 
equation of the components of density matrix is rewritten as 

dtpa = {Kp + Kl3)Pp- (19) 

The matrix h'^^j must obey some restrictions. The probability is conserved, entropy 
must be real and never decrease, which imply that h'^Q = h'^f^ = and /i^^ is negative 
semidefinite. The energy conservation in neutrino oscillation demands that h'^Q = H'q^ = 
have to be imposed. The /i^^ is thus parameterized as 
















—a 


-p 








-P 


-7 





V 








J 



(20) 



a > 0, 7 > 0, a7 > /?^. 
By solving the Eq.(19) we get the density matrix 

~ \ lsin29{pi + ip2) sin^e J ' ^ ^ 
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where the pi and p2 have different forms due to the magnitude relation between the 
difference of neutrino mass square and the QMV parameters, which will be discussed as 
follows. 



(i) If the condition 



{a - -if + 4/32 < 



A 

'e 



is satisfied, then p has a oscillation-like solution 

Pi = cosSit e 2^*, 

1 — 7 



a+j_. 



P2- A 



^-3' 2 



cos6it — 6isinSit)e 2 , 



where 



c _ A 



1 - 



{a - 7)2 + 4/32 



The oscillation probability is 



P{j^e i^/i) = ]^sin^2e{l - cosSit ■ e""^*). 



(22) 



(23) 



(24) 



(25) 



(ii). If the relation 



(a -7)2 + 4/32 > 



A 

'e 



is valid, the p has an exponential-like solution 



Pi = [exp i^t) + exp (-^^)]e 2 



(26) 



P2 = A 



-^^[(a - 7 + 52)exp(^t) + (a - 7 - 52)exp(-^i)]e 



here 



V -7)2 + 4/32 -(|)2. 



The probability of neutrino oscillation becomes 

^(^e ^ ^^) = -sm22^[l - -(exp( L U) + exp{ L U))]. 



(27) 



(28) 



(29) 



It is easy to show that Trp^ ^ the pure state may evolve into mixed state. With 
the q;,/9, 7 — > 0, the p{t) return to Eq.(6) in these two situations, quantum mechanics 
recovers. 
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4 QMV for neutrino oscillation in matter 



In the situation of neutrino oscillation in matter, we add the parameterized QMV 
term of Eq.(20) into the Hamiltonin of Eq.(15) and get the evolution equation of the 
components of density matrix 



dt 



( Po\ 




/ 


Pi 




-a 






-A- 


\ P3 ) 




^ 



\ 

-13 

-7 

/ 



( Po \ 
Pi 
Pi 
\pz J 



(30) 



Taking the adiabatic approximation and the initial condition that only electron neutrinos 
exist in the system into account , solving differential equations Eq.(30), we obtain the 
density matrix 

/ s _ / cos^O \sin29{pi - ip2) 



sin^O 



(31) 



where 9 is the effective mixing angle of neutrinos in matter at the initial moment. The 
Pi and p2 are given as below. 



(i) . In the situation of 



then 



P2 



\-(3' 2 



{a - -if + 4/3^ < 4A^ 

Pi — cosi / 52,dt') • 2 *^ 
Jo 

{f 5^dt')-5^sin{ f's^dt')] 
Jo Jo 



(32) 



1 .a — '-f 
;[ — — COS 



e 2 



where 



53 = -V4A2-(a-7)2-4/32. 



The neutrino oscillation probability is calculated as 

Pii^e ^ i^u) = ^(1 - COS 2^ cos 2^) - \sin2esin29 ■ cos( f 6^dt')e~^ 
2 2 Jo 



(ii) In another situation 
then 



P2 



where 



2(A-/3) 



(a - -if + 4/3^ > 4A^ 

Pi = -coshi I 5idt')e ~* , 
2 Jo 

q; — 7 — 5Aexpi— / 5idt')]e~ 
Jo 

6, = ^^(a -7)2 + 4/32 -4A2. 



(33) 
(34) 

(35) 
(36) 



1 /•* 

[(a — 7 + 64)exp( / 64dt') + 
Jo 



Q+7 , 



(37) 
(38) 
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The neutrino oscillation probability becomes 



Pii^e z^m) = ^(1 - cos2ecos2e) - ^sin2esin2ecosh{J^ 5idt')e-^K (39) 

It is not difficult to see that with the QVM parameter a, /3, 7 0, the probability 
Piye — ^ J-'fi) back to Eq.(18) in the two situations, which is adiabatic MSW solution 



o2 



of solar neutrino problem. It is easy to calculate that Trp < 1. In this case, pure state 
may evolve into mixed state, and thus quantum mechanics is violated. 



5 An estimation for the bounds on QMV parame- 
ters 

We have derived the formulas of QMV effect in neutrino oscillation, which is to be tested 
in experiments. A complete analysis requires a detailed understanding of all neutrino 
experiments which goes beyond the scope of this paper. We present here only a rough 
estimate of maximum order of magnitude of the QMV parameters as an illuminating 
example of testing quantum mechanics in neutrino oscillation. It is estimated theoreti- 
cally that the maximum possible order of magnitude for QMV parameters a , , or 7 
is Q(-E^/mp;)|]10||, which in K^K^ system is ~ lO'^^Gel^, where £" is a typical energy 



scale in the system under discussion, and Mpi is the Plank energy scale. The bounds 
on these parameters from K^K^ experiments have been obtained in Refs.[^,§] and [ |T0[] , 
the last one gives 

ctKORO < 4 X 10"^^ GeV, \/3koko \ < 3 x 10"^^ GeV, -f^oKO < 7 x lO^^^ GeV. (40) 

Similary, in neutrino system, it is expected theoretically that the order of magnitude of 
the maximum parameter is 0{El/mpi), where E^, is a typical energy scale of neutrino 
system. This is the order of ~ 10~^^GeV for solar neutrinos or reactor neutrinos. 
In general, we may reasonably assume that the order of the maximum parameters in 

neutrino system is 0{ ^pT" "^ g^o-o ~ ^ e^o-o Y'^k'^Ro, and that the magnitude relation 
among these parameters is retained in neutrino system, i.e., a ^ \P\ ^ 7. 

In the formulas of neutrino oscillation probability, the QMV parameters appear al- 
ways in company with time t, for example at. There is no enough large time to have the 
at ~ 1 in current territorial experiments. Therefore, we discuss here only solar neutrino 
experiments because of large neutrino propagation time. For the situation of neutrino 



oscillation in vacuum which is a possible solution of solar neutrino problem[|T6|, From 
Eq.(25), we get 

^ 2, .sm^2e-2P. 

where (and hereafter) the P(t'e J^^t) is simply written as P. If the condition Eq.(26) 
is satisfied, we get 

« + 7>§- (42) 
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/^From above arguments of the order of magnitude of the maximum parameters and 
the order of maximum parameter of K^K^ system in Eq.(40), we get an upper hmit of 
another allowed region 

(a + 7)max<1.6xlO-20G'eV^ (43) 

Two allowed regions of the a + 7 are decided by Eq.(41)-(43). 

We studied the upper limit of a + 7 given by Eq.(41) vary with the sin^ 2^, the 
numerical results is shown in Fig.l. The lower and upper limits of the second allowed 
region are also shown in same figure. The neutrino oscillation probability P is taken 
from the Kamiokande0 experimental data and the prediction of standard solar model 
(SSM)|T^. The region of sin'^26 is given by recent fit]l8| of world solar neutrino exper- 
iments. The A is taken as the center value of A ~ 6 x lO^^^GeV^, and E ^ lOMeV is 
used. 

In other situation, we discuss neutrino oscillation in matter, which is another possible 
solution (MSW effect) ||15| of solar neutrino problem. /^From Eq.(35) we get 

2, A - cos29cos29 - 2P , ,^ , 

a + 7< — In = . (44) 

' - t ' sin2esin2e ' ^ ^ 

/^From Eq.(37) and Eq.(21), we get 

a + 7 > 2A. (45) 
Two allowed regions of a + 7 are given by Eqs.(43)-(45) for the case of including matter 



effect . iFrom the Kamiokande experimental data, SSM predictions, as well as the fit|118 
of experimental data to the adiabatic MSW solution of solar neutrino problem, we obtain 
the numerical result of variation of the parameter a + 7 with the sin^29, which is shown 
in Fig. 2. In numerical calculation, we used the center value of the A given in Ref.[0, 
A ~ 1.6 X 10~^el^^, and E ~ lOMeV. We also apply the average value of the sin^O 
over the radius of the Sun. The electron distribution in the Sun is taken from Ref . [pPj] . 

It should be emphasized that our numerical result is only an example far from an 
exact analysis. The aim is to show how the information of bounds on QMV parameters 
is extracted from neutrino oscillation experiments. The value of the bounds are different 
if the different experimental data are used. A complete numerical analysis based on our 
analytical formulas is necessary. It is noted that the QMV for neutrino system in vacuum 



was discussed in Ref. [19]. However, there is no analytical result and no any information 



about bounds on parameters in their disscussions. Moreover, what discussion in Ref. [ff^ 



is the case of neutrino oscillation in vacuum, but the neutrino mixing parameters used 



in their calculation are taken form the values of MSW effect given in Ref. |^ , there is 
no any room of small mixing angle to be allowed by experiments for vacuum oscillation 
solution of solar neutrino problem. 



6 Conclusions 

Neutrinos may be a useful system to precisely test quantum mechanics, the Liouville 
equation that describes the time evolution of density matrix of neutrino system can be 
modified by adding an extra term, which may be parameterized as a 4 x 4 matrix. The 
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modified evolution equation of density matrix can be solved analytically if neutrinos 
propagate in vacuum or if neutrinos propagate in non-uniform matter but adiabatic 
condition is available. The analytical expressions of neutrino oscillation probability with 
or without QMV effect have been derived. Based on theoretical analysis and neutrino 
experiment results we have extracted the bounds on the QMV parameters. Two allowed 
region of a + 7 have been obtained. It is expected that more precise restrictions on QMV 
parameters may be obtained from future neutrino experiments and complete numerical 
analysis. 
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Appendix A: Generalized adiabatic condition 

In this appendix we will demonstrate that if the adiabatic condition is valid the term 
being proportional to {dX/dt) can be neglected when we solve the evolution equation of 
the components of density matrix . /^From the Eq(14) we get 

IdX (A/A) - cos 2^ 1 dX 



Xdt . /^cos 29 - A/A f + sin^ 29 ' ^X dt 



(Al) 



We define two dimensionless parameters C and G as 



1 dX , , , 

where E is a typical enegy scale in the system under discussion, which is taken as 
neutrino energy in the system considered here. 

(.3) 

where m^i_2 are the eigenvalues of the matrix defined in Eq.(lO). The adiabatic 
condition is expressed |T^ as G <^ 1. From Eqs.(ll) and (Al)-A3), we get 

C _ A I sin 2^ cos 2^1 

G~E^ sm^29 ' ^ ^ 

For the problems related to current neutrino experiments, we estimate the maximum 
valueoftherighthandedof Eq.(A4). Amax < lOeV^ , Emin ~ IMeV , and (| sin2ecos2^|/sin2 2^), 
10^ , We thus get 

^ < 10-^ (A5) 

If adiabatic condition is valid, G <^ 1 , we have C ^ 1 . Therefore, we get the generalized 
adiabatic conditions in density matrix description of neutrino system, we can safely 
neglected the term containing (dX/dt) during we solve the differential equation of density 
matrix as long as the adiabatic condition is satisfied. Under adiabatic approximation 
the analytical solution of the evolution equation can be obtained. 
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Figure 1: 
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Figure 2: 
Figure Caption 

Fig.l The limits of the parameters a + 7 vary with the sin^2d for neutrino oscillation 
in vacuum. The curve III is an upper limit for the case of (a — 7)^ + 4/9^ < ^. The line 
II is a lower limit for the case that above condition is not valid. The line 1 is an upper 
limit from the theoretical arguements and the experimental upper limit of the maximum 
parameter in K^K^ system. The arrows indicate the two allowed regions. 
Fig. 2 Same as Fig.l except for the neutrino oscillation in matter. 
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